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We report on the first observations of B° Di (2430)°u;, B° —>■ Di(2420)°u> and B° —> 
D|(2460)°cj decays. The B° -> D* + p(1450)~ decay is also observed. The branching fraction 
measurements are based on (771.6 ± 10.6) x 10 6 BB events collected at the T(4S) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e“ collider. The fractions of longi¬ 
tudinal polarization of the D** states as well as partial wave fractions of the Di(2430)° are ob¬ 
tained. We also set a 90% confidence level upper limit for the product of branching fractions of 
B(B° —> D* + 6i(1235)“) x £>(&i(1235) _ —> uott~). The measurements show evidence of nontrivial 
final-state interaction phases for the p-meson-like amplitudes. 
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I. INTRODUCTION 


Orbitally excited states of the £> meson (£>** 
states) provide a good opportunity to test heavy 
quark effective theory (HQET) [l] and QCD sum 
rule predictions [2]. The simplest system consists of 
a charm quark and a light antiquark in an orbital 
angular momentum L = 1 (P-wave) state. Four 
such states are expected with spin-parities J p = 0 + 
(j = 1/2), 1+ (j = 1/2), 1+ ( j = 3/2) and 2+ 
(j = 3/2), where j is the sum of the light quark 
spin and angular momentum £. All these states have 
been discovered [||. They are Dp (2400), £>i(2430), 
£>i(2420) and £>2(2460). The conservation of par¬ 
ity and angular momentum in strong interactions 
imposes constraints on the decays of £>** states to 
£>(*) 7r. The j = 1/2 states are predicted to decay 
mainly through an S- wave: £>q(2400) —> Dn and 
£>i(2430) —> D* 7T. The j = 3/2 states are expected 
to decay mainly through a £>-wave: £>i(2420) —> 
D*n and £>2(2460) — > Dn and D*n. The j = 1/2 
states with £ = 1 are expected and proven to be 
broad (hundreds of MeV/c 2 ), while the j = 3/2 
states are expected and proven to be narrow (tens 
of MeV/c 2 ). The BABAR [|| and LHCb [bj collab¬ 
orations have discovered other excited £> mesons in¬ 
terpreted as n£ = 2 S and n£ = 1£> states as well as 
a possible superposition of several nL = IP states, 
where n is the radial quantum number. 

Since HQET is violated, the physical £>** state 
with J p = 1 + can contain admixtures of the states 
with j = 1/2 and j = 3/2 Q. 

A similar spectroscopy exists for the D s j states 
0. However, the observed masses for the £>* 0 (2317) 
and £> s i(2460) resonances with j = 1/2 are signifi¬ 
cantly smaller than predicted [jj. The £> s i(2536) and 
££i(2460) states with J p = 1 + can mix with each 
other. This effect is observed in an angular analysis 
of the £> s i(2536) + —» D* + Kg decay Q. 

Precise knowledge of the properties of the £>** 
states is important to reduce uncertainties in the 
measurements of the semileptonic decays and thus 
in the determination of the Cabibbo-Kobayashi- 
Maskawa matrix elements \V c b\ and \V u b\ 

The £>** mesons have been observed in both 
semileptonic [10] and hadronic B decays [lTI - fl^ |. 
The recent LHCb study [l4] shows the first obser¬ 
vation of the B° —> £>3(2760) + 7r _ decay as well. 
The dynamic properties of £>** production are de¬ 
termined by the Wilson operator product expansion 
m . In color-favored B° —> D** + n~ decays [111 
dominance of the narrow £>** states over the broad 
ones is observed. A study of B~ —> D**°n~ decays 
IT3 with the color-favored and color-suppressed pos¬ 
sibilities shows approximately equal production of 


the broad and narrow £>** mesons. It can be ex¬ 
plained by a significant suppression of the narrow 
states in the color-suppressed channel. Calculations 
based on HQET and quark models [l6| predict such 
suppression. 

In this paper, we perform an amplitude analysis of 
the B° —► D* + uin~ decay to measure the decay frac¬ 
tions to £>** states produced via the color-suppressed 
channel [Fig. [Ha)] and to study the £>** properties. 


This decay is sensitive not only to the vector 
£>i(2430)° and £>i(2420)° states but also to the ten¬ 
sor £>2(2460)° state. Although the B° decay to 
£>2(2460)°w is prohibited under the naive factor¬ 
ization hypothesis, it can nevertheless be produced 
via final-state interactions (FSI) and/or nonfactoriz- 
able contributions. In soft-collinear effective theory 
(SCET), color-suppressed decays to the £>|(2460)° 
state can receive a factorizable contribution at the 
leading order in A /ms G2- This mechanism leads 
to the equality of branching fractions and strong 
phases in the decays B° — > £>2(2460)°M and B° —> 
£>i(2420)°M, where M is a light meson. Possible 
deviations from this equality can be attributed to 
subleading effects 0- A discussion of £>** produc¬ 
tion in hadronic B decays can be found in Ref. [l8| . 

The color-favored mode of the studied decay 
[Fig. mb)} i s saturated by light con resonances. 
Hadronic weak currents can be classified as either 
first- or second-class, depending on the combination 
of spin J and the P- and G-parities of the cun sys¬ 
tem [l9[ . In the standard model, first-class currents 
(FCC) have J PG = 0 ++ ,0 , 1 H or 1 and are 

expected to dominate. Second-class currents (SCC) 
have J PG = 0 +_ ,0“ + ,l ++ or 1 and are associ¬ 
ated with a decay constant proportional to the mass 
difference between the up and down quarks. Thus, 
they are expected to vanish in the limit of perfect 
isospin symmetry. The decay B° —> D* + uin~ is ex¬ 
pected to proceed predominantly through the FCC, 
mediated by p-meson-like resonances, such as off- 
shell p(770)“ and p(1450) - . In contrast, the SCC 
may be mediated by &i(1235)“. SCC searches have 
been performed extensively in nuclear /3 decays [2C} 
and r decays [2l|, with no evidence found. 

The structure of the p-meson-like states is not yet 
completely clear. The p(1450) has a mass consistent 
with that of a radial 2 S excitation [22] but its decays 
show characteristics of hybrids [23[ and suggest that 
this state may be a 2S'-hybrid mixture 24]. The 
observation of the p(1450) in B -meson decays and 
the study of its interference with the p(770) would 
lead to a better understanding of the properties of 
the p-meson-like states. 

Another aim of this study is a test of the factoriza¬ 
tion hypothesis in the £>** production region. The 




Figure 1: (color online), (a) Color-suppressed and (b) color-favored tree diagrams for the production of D** and uin 
states in B° —» D* + uon~ decays, respectively. 


factorization hypothesis, widely used in heavy-quark 
physics for hadronic two-body decays, assumes that 
two hadronic currents may be treated independently 
of each other, neglecting FSI. The factorization can 
be tested by examining the polarization in H-meson 
decays into two vector mesons. The idea is that, 
under the factorization, certain hadronic decays are 
analogous to similar semileptonic decays evaluated 
at a fixed value of the momentum transfer, q 2 = 

[25l | . Based on the polarization measurements of the 
decays B° — > D*°uj [26} and B — > (j>K* [ 13 ], we can 
conclude that nonfactorizable QCD effects are es¬ 
sential in color-suppressed decays. The significant 
transverse polarizations measured in these decays 
may arise from the existence of effects from non¬ 
trivial long distance contributions, as predicted by 
SCET studies [28} . The longitudinal polarizations of 
similar decays B° -> £>i(2430)V B° Di(2420)°w 
and B° —> £>2(2460)°w are measured in our study. 

The studied decay has been first observed by the 
CLEO [13| and BABAR. f3fj collaborations, the lat¬ 
ter finding an enhancement in the D*tt mass broadly 
distributed around 2.5 GeV/c 2 . 

II. EXPERIMENT AND DETECTOR 

This study uses a data sample containing 771.6 ± 
10.6 million BB events collected at the T(45) 
resonance with the Belle detector at the KEKB 
asymmetric-energy e + e _ collider j3lj. The Belle de¬ 
tector, which is a large-solid-angle magnetic spec¬ 
trometer based on a 1.5 T superconducting solenoid 
magnet, consists of several subdetectors. 

Charged particle tracking is provided by a 4-layer 
silicon vertex detector (SVD) and a 50-layer central 
drift chamber (CDC). The charged particle accep¬ 
tance covers laboratory polar angles between 6 = 
17° and 150°, corresponding to about 92% of the 
total solid angle in the e + e _ center-of-mass (c.m.) 
frame. 

Charged hadron identification is provided by the 


ionization energy-loss cLE/dx measurements in the 
CDC, an array of aerogel threshold Cherenkov coun¬ 
ters (ACC), and a barrel-like arrangement of tirne- 
of-flight scintillation counters (TOF). The informa¬ 
tion from these three subdetectors is combined to 
form likelihood ratios (PID), which are then used 
for pion, kaon and proton discrimination. An elec¬ 
tromagnetic calorimeter (ECL), comprised of 8736 
CsI(Tl) crystals and covering the same solid angle as 
the charged particle tracking system, serves for the 
detection of electrons and photons. Electron iden¬ 
tification is based on a combination of dE/dx mea¬ 
surements in the CDC, the response of the ACC and 
energy-to-momentum ratio of an ECL shower with 
a track as well as a transverse shape of this shower. 
An iron flux-return located outside of the coil (KLM) 
is instrumented to detect mesons and to iden¬ 
tify muons. The detector is described in detail else¬ 
where [H}. 

The EvtGen event generator, [33] with PHOTOS 
{Hj for radiative corrections and a GEANT-based 
Monte Carlo (MC) simulation [35} to model the re¬ 
sponse of the detector and determine the accep¬ 
tance, are used in this analysis. The MC simu¬ 
lation includes run-dependent detector performance 
and background conditions. 

III. EVENT SELECTION 

Candidate B° —> D* + ojtt~ events as well as 
charge-conjugate combinations are selected. The 
D* + candidates are reconstructed in the D* + —> 
D°n + mode. The D° candidates are selected us¬ 
ing the D° —> K~tt + mode. Other D° decay modes, 
which lead to significantly smaller signal-to-noise ra¬ 
tios, are not used in this analysis. The to candidates 
are reconstructed in the to —> ir + ir~ir° mode. 

Charged tracks are selected with a set of track 
quality requirements based on the average hit resid¬ 
uals and impact parameters to the interaction point. 
To reduce the low momentum combinatorial back- 
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ground, we also require that the track momentum 
transverse to the beam direction be greater than 100 
MeV/c for all tracks except for the slow pion candi¬ 
date in the D* + —> D°tt + decay, for which we apply 
a looser cut of 50 MeV/c. 

A PID requirement is applied for kaon candidates 
but not for pion candidates. The kaon identification 
efficiency is about 90% and the pion misidentffication 
rate is less than 10%. All tracks that are positively 
identified as electrons are rejected. 

Photons are identified as ECL clusters that are not 
associated with charged tracks and have a minimum 
energy of 70 MeV in both the barrel and endcap 
regions. 

D° candidates are reconstructed from com¬ 

binations with an invariant mass within 15 MeV/c 2 
of the nominal D° mass [3]. This window corre¬ 
sponds to approximately ±3 times the mass resolu¬ 
tion. D* + candidates are selected by combining D° 
candidates with an additional track, assumed to be a 
7r + . The mass difference — m^o is required to be 
within 2 MeV/c 2 of its nominal value; the resolution 
of this quantity is about 0.5 MeV/c 2 . 

Neutral pion candidates are formed from pho¬ 
ton pairs that have an invariant mass within 11.25 
MeV/c 2 of the nominal tt° mass, which corresponds 
to about ±2.5 times the reconstructed mass resolu¬ 
tion. To reduce the combinatorial background, the 
total energy of the photons is required to be greater 
than 250 MeV. 

The u> candidates are formed from a pair of 
oppositely-charged tracks, assumed to be a 7r + 7r _ 
pair, and a n°. The invariant mass of the 7r + 7r _ 7r° 
combinations is required to be within 73.5 MeV/c 2 
of the nominal w mass. This very loose cut retains 
sideband candidates for background estimation. The 
instrumental resolution on the ui candidates is about 
7.3 MeV/c 2 . 

To reduce the number of false to candidates formed 
from random combinations of pions, we impose an 
additional requirement in the to Dalitz plane, moti¬ 
vated by the oj decay dynamics and spin-parity in 
the B° —> D* + lutt~ decay j36|. We define two or¬ 
thogonal coordinates X = 3 T 0 /Q — 1 and Y = 
V3(T+ - T-)/Q, where T±. o are the kinetic energies 
of the pions in the u> rest frame and Q = Tb±T_ ±7+ 
is the energy release in the u> decay. Further we de¬ 
fine a variable r properly scaled to the kinematic 
limit as 

r = — ’ (1) 

where r b is the distance from (0, 0) to the boundary 
in the direction of (A', V). Since the Dalitz plot den¬ 
sity peaks at r = 0 for the ui signal, we impose the 


requirement r < 0.75. This requirement eliminates 
about 41% of the background while retaining about 
84% of the signal. In Fig. [2] we show the simulated 
(X, Y) Dalitz plane of the uj signal events and the 
restriction on r variable. 
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Figure 2: (color online). Simulated (X, Y) Dalitz distri¬ 
bution of the w signal events. The curve bounds the area 
selected for further study. 
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B candidates are reconstructed by combining 
a D* + candidate, an ui candidate, and an addi¬ 
tional negatively charged track. All B candidates 
are identified using two kinematic variables: the 
energy difference A E = JV \/|p*| 2 c 2 + m?c 4 — 
fi^eam and the beam-constrained mass M^ c = 
t/^beam/ c4 ~ I Si P*| 2 /c 2 , where the summation is 
over all particles forming the B candidate, p* and 
mi are their three-momenta and masses, respec¬ 
tively, and £'b eam is the beam energy. All quanti¬ 
ties are defined in the e + e _ c.m. frame. We select 
events with a tight cut on Mbc of 5.2725 GeV/c 2 < 
Mbc < 5.2845 GeV/c 2 , corresponding to about ±2 
times the mass resolution, and a loose cut on A E 
of |AF?| < 0.22 GeV. To suppress possible contin¬ 
uum events (e + e _ — > qq, where q = u,d 7 s,c ), we 
limit the angle between the thrust of the B candi¬ 
date and that of the rest of the event by requiring 
| COS ©thrust | <0.8 H3- 

In this study, we perform an amplitude analysis 
that accounts for the kinematic properties of the de¬ 
cay matrix element. The matrix element should be 
symmetrized relative to the exchange of two identi¬ 
cal particles in the final state (two tt~ mesons in our 
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decay mode fUJ) according to the identity principle. 
Such symmetrization leads to an interference term 
in the squared matrix element. This term consists of 
two lo decay amplitudes with different 7 r + 7 r _ 7 r° com¬ 
binations in the D* + n + n~n° tt~ final state. Since 
the lo is a relatively narrow resonance, the interfer¬ 
ence term is essential only in the overlapping region 
of the 7 T + 7 r - 7 r° invariant masses, which is of the or¬ 
der of the lo width. To correctly describe the angular 
distributions in this interference region, the internal 
degrees of freedom of the lo decay should be taken 
into account 136]. However, due to lack of statistics, 
which prevents a full analysis in such a case, we ex¬ 
clude this interference region without significant loss 
of statistical power. 

In order to reduce smearing from detector resolu¬ 
tion, a simultaneous fit constraining the 77 , K~ 7 r + , 
D°n + and D*+ lott invariant masses to match the 
known 7 r°, D °, D* + and B° masses, respectively, is 
performed. The 7 r + 7 r _ 7 r° invariant mass is not con¬ 
strained to the lo mass in the fit because of the non- 
negligible width of the lo meson. 

There are events for which two or more candidates 
pass all the selection criteria. According to MC simu¬ 
lation, this occurs primarily because of the misrecon- 
struction of one of the pions from the lo —> 
decay. To ensure that no B decay is counted more 
than once, a best-candidate selection is performed 
based on a y 2 defined as the sum of three terms. 
The first determines the deviation of the 7 r° invariant 
mass from its nominal value, the second represents 
the deviation of Mbc from the nominal B° mass and 
the third uses the distribution of the difference be¬ 
tween the ^ coordinate at the interaction point of the 
track corresponding to the primary pion ( 7 r~) from 
the B° signal decay and the average z coordinate for 
the tracks corresponding to the decay products (K~ 
and 7 r + ) from the D° meson decay. We retain only 
the z coordinate information because B mesons are 
boosted along z and the vertex resolution is worse 
in that direction. We omit the lo candidate mass in 
this procedure in order to avoid any bias in the u> 
mass distribution since this distribution is used ex¬ 
tensively for the background description. 

The signal sample is composed of two 
components—correctly reconstructed (CR) and 
self cross-feed (SCF)—that are distinguished by 
whether or not the kinematic variables of the 
D* + lott~ decay are well reconstructed. MC simula¬ 
tion shows that the SCF component predominantly 
occurs due to the combinatorial background for 
the lo. To define the CR and SCF components, we 
use the following y 2 describing the deviation of the 
reconstructed momenta of the final particle system 


(rec) from the generated momenta (gen): 

3 /•„(*) _ r (0 ^2 

2 _ V"' V^fcgen J 'fcrec/ 

X ~ 2 -^ 2 -^ 2^(0 ) ’ ' ' 
i k= 1 ® v*^fcgen/ 

where x^ \ 3 = {P^\ are the spherical mo¬ 

mentum coordinates of the ith particle in the final 
state, cf(x^) is the corresponding detector resolu¬ 
tion, and the summation is over all tracks and 7r° 
forming the B candidate. We choose to define the 
CR (SCF) component by the condition y 2 < C 
(y 2 > C). The value of C — 300 is determined 
by examining the shapes of the distributions of the 
difference between the reconstructed and generated 
kinematic variables. Variations of the value of C are 
considered as a source of systematic uncertainty. 

Figure [3] shows the distribution of the selected 
events in the (AE, M(tt + ir~n°)) plane, where we de¬ 
fine the following four regions to distinguish between 
signal and background: 
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Figure 3: Distribution of A E versus M('k + tv~'k°) for the 
selected B° —► D* + ion~ candidates. The signal Region 
(I) and sideband Regions (II, III and IV) are shown. A 
clear correlation between the variables is seen in Region 

I. 


I -> \AE\ < 34 MeV; 

|M( 7 r + 7 r _ 7 r°) — m^l < 21.25MeV/c 2 , 
II -> 66 MeV < \AE\ < 198 MeV; 

|M( 7 r + 7 r _ 7 r°) — m u \ < 21.25 MeV/c 2 , 
III -> |A£7| < 34MeV; 
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\M(tt + tt 7t°) — m w \ € [34; 68] MeV/c 2 , 

IV -> 66 MeV < \ AE\ < 198 MeV; 

|M(7r + 7r _ 7r°) — m^\ £ [34; 68] MeV/c 2 . 

Here m u is the nominal u> mass. Region I is the 
signal region while the others are sideband regions. 
A clear correlation between the A E and M (7r + 7r _ 7r°) 
variables is seen in Region I due to the experimental 
resolution. The signal window for the uj invariant 
mass corresponds to ±2.5 times the world-average u> 
width of 8.5 MeV/c 2 . 

Figure S] shows the AI(n + TT~n°) distributions in 
the A E signal and sideband regions defined above. 
The curve corresponds to the sum of a Voigtian func¬ 
tion (the convolution of a Breit-Wigner function with 
a Gaussian function) and a linear background func¬ 
tion. The uj mass, the Gaussian resolution a , and the 
parameters of the linear function are free in the fit 
but the Breit-Wigner width is fixed to the world- 
average decay width of the uj @. The difference 



M(7iVti 0 ) (GeV/c 2 ) 

Figure 4: M(n~Tr + n°) distribution of the B° D* + um~ 
candidates in the A E signal region (points with error 
bars) and sideband (hatched histogram). The sideband 
distribution is normalized to the size of the A E signal 
region. The M(n + n~n°) signal region and sideband are 
indicated by the vertical lines. The curve is the result of 
the fit described in the text. 

between the number of observed events away from 
the M(tt + tt~tt 0 ) peak and the number of events pre¬ 
dicted from the A E sideband is explained by the 
A E peaking background component, corresponding 
to B° —> D* + TT + n~Tr 0 Tr~ decays. 


IV. TOTAL BRANCHING FRACTION 

The signal yield is obtained from a binned y 2 
fit to the A E distribution using a function describ¬ 
ing the CR and SCF components together with a 
smooth combinatorial background. Since the B° —> 
D* + 7r + 7T _ 7T° n~ events observed in Fig. 3 produce 
a peak in A E, the fit is performed separately in 
the M{tt + tt~tt°) signal and sideband regions defined 
above on the (AE, M(n + Tr~TT 0 )) plane. MC simula¬ 
tion shows that these events have the same shape 
as the CR component. In the fit, the CR com¬ 
ponent is described by a double-Gaussian function 
with distinct means and widths, the SCF compo¬ 
nent is described by the sum of a Gaussian function 
and a second-order polynomial, and the combinato¬ 
rial background is described by another second-order 
polynomial. The means, widths and relative normal¬ 
izations of the CR and SCF functions are fixed to the 
values obtained from the signal MC simulation, while 
the signal normalization and the parameters of the 
polynomial background function are treated as free 
parameters. The differences between MC and data 
values for the fixed parameters in the fit are found to 
be within MC statistical errors. The fit results are 
shown in Fig. [5] in both the r°) signal and 

sideband regions. The fitted signal yield is found to 
be 919 ± 37 for the M(7r + 7r _ 7r°) signal region and 
157 ± 21 for the sideband region. The final yield 
Ns = 821 ± 39 is computed as the difference be¬ 
tween these two yields, taking into account the ratio 
of 5/8 between the widths of the M(7r + 7r _ 7r°) signal 
and sideband regions. 

The fraction of neutral B mesons decaying to the 
studied final state is expressed as 

B = 131 

es^VB^ec’ 1 ’ 

where es = (2.11 ± 0.02)% is the detection efficiency 
determined from a MC simulation that uses a Dalitz 
plot distribution generated according to the signal 
model described below, r/ = 0.941 ± 0.029 is the ef¬ 
ficiency correction factor that accounts for the dif¬ 
ference between data and MC and obtained from 
the momentum-dependent corrections for the n° and 
slow pion from the D* decay and the PID corrections 
for the kaon, Nb = (771.6 ± 10.6) x 10 6 is the to¬ 
tal number of neutral B mesons in the data [39| and 
£> sec = (2.32±0.04)% is the product of the secondary 
branching fractions. Using Eq. ©, we obtain 

B = (2.31 ± 0.11 (stat.) ± 0.14 (syst.)) x 10“ 3 , 

which is consistent with the CLEO value (29j within 
1.2(7 and the BABAR value f30] within 1.5cr. The 






















AE (GeV) 


AE (GeV) 


Figure 5: AE distributions of the B° —¥ D* + con~ candidates in the (a) signal and (b) sideband regions of M(tv + tt~ty 0 ). 
The hatched histogram in (a) represents the M(tv + tv~tv°) sideband normalized to the size of the signal region. The 
AE signal region and sideband are indicated by the vertical lines. The curves are the results of the fit described in 
the text. 


Table I: Sources of relative systematic error in the 


branching fraction measurement. 

Source Error (%) 

Signal yield, Ns 

— M(n + n~n°) signal region 1.3 

—Definition of SCF and CR components 0.9 
—AE signal shape 2.2 

—AE background shape 1.3 

Signal efficiency, es 

—Track reconstruction efficiency 3.9 

—7T° reconstruction efficiency 2.3 

—Kaon identification efficiency 0.9 

— B° signal decay model 1.1 

—MC statistics 0.8 

Number of neutral B mesons, Nb 1.4 

Secondary branching fractions, B sec 1.7 

Quadratic sum 6.1 


total systematic error of 6.1% summarized in Table U 
arises from the following sources: 

(i) An uncertainty of 1.3% due to the choice of 
the signal window for the M(tt + ir~n°) invari¬ 
ant mass is estimated by reducing the size of 
the window from 21.25 to 12.75 MeV/c 2 . The 
reduced window corresponds to 1.5 times the 
world average w width. 


(ii) An uncertainty of 0.9% related to the definition 
of the SCF and CR components is estimated by 
changing the requirement on the y 2 defined in 
Eq. © to C = 200 or C = 400. 

(iii) An uncertainty of 2.2% related to the A E shape 
description is estimated by varying the shape 
parameters fixed from MC simulation in accor¬ 
dance with their MC statistical errors. 

(iv) An uncertainty of 1.3% due to the background 
description in the AE shape is estimated by 
adding higher-order polynomial terms or keep¬ 
ing a linear term only. 

(v) A dominant uncertainty of 3.9% is assigned to 
the total reconstruction efficiency of all charged 
tracks in the decay. For a single track, this un¬ 
certainty depends_on the transverse momentum 
Pt of the track [4p| . For low momentum tracks 
(with pt < 200 MeV/c), it is estimated using 
the decays B° —> D*~tt + and B + —> D*°n + ; 
for high momentum tracks, a study of the track¬ 
ing efficiency is based on partially reconstructed 
D*+ —> D°(KgTT + n~)TT + decays. The total 
tracking error is the linear sum of the errors 
corresponding to the individual tracks. 

(vi) An uncertainty of 2.3% in the reconstruction 
efficiency of neutral pions is estimated using the 
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r —» 7r 7 t°z/ t branching fraction and events 
where the other r decay is tagged [iH] . 

(vii) An uncertainty of 0.9% in the efficiency of 
the kaon particle identification requirement is 
obtained using a control sample of D* + —> 
D 0 (K~n + )Tr + decays [42j . 

(viii) An uncertainty of 1.1% is assigned due to the 
model dependence of the signal reconstruction 
efficiency. The signal model with the best de¬ 
scription of the data is constructed in Secs. IV.51 
and IV.Cl in the frame of the amplitude analysis. 
The model parameters obtained from the fit in 
Sec. IV.Cl have statistical uncertainties. These 
are propagated as a systematic uncertainty on 
the signal efficiency, taking into account the full 
covariance matrix. 

(ix) A binomial uncertainty of 0.8% due to the lim¬ 
ited Monte Carlo sample size arises in the effi¬ 
ciency calculation. 

(x) An uncertainty of 1.4% in the number of B 
mesons is estimated from Ref. [43| . 

(xi) An uncertainty of 1.7% is associated with the 
measured branching fractions of the D*, D and 
w 0. 


V. AMPLITUDE ANALYSIS 

To study the resonant structure of the B° —> 
D* + ujit~ decay, we perform an amplitude analysis. 
Using an unbinned likelihood method, we simultane¬ 
ously fit the data in the six-dimensional phase space 
according to Ref. [3(|. 

We define two sets of kinematic variables: 
[M 2 (uott), cos 0i, </>i, cos/3i, and cos£i] and 
[M 2 (D*tt), cos 9-2 , 02, cos/32, 02 and cos^], cor¬ 
responding to the ujtt and D** productions (color- 
favored and -suppressed diagrams of Fig. [[]), respec¬ 
tively. 

The masses M{ujtt) and M(D*tt) are the invari¬ 
ant masses of the ujtt and D*tt combinations. The 
angular variables, [cosfi 1 !, 0i, cos/3i, 0i and cos£i], 
describing ujtt production, are defined in Fig. [G] The 
polar and azimuthal angles, 9\ and 0i, defined in the 
uj rest frame, are the angle between the normal to 
the lo decay plane and the ujtt direction, and the an¬ 
gle between the R-decay plane and the plane formed 
by the n u and ujtt directions, respectively. The polar 
and azimuthal angles, /3i and 0i, defined in the D* 
rest frame, are the angle between the D and the ujtt 
flight directions, and the angle between the B- and 
D*- decay planes, respectively. The polar angle £i is 


the angle between the D* and to flight directions in 
the ujtt rest frame. 

The angular variables, 02 and 02 as well as /?2 and 
02, describing the D** production, are defined in 
the same manner as angles for the ujtt production 
but with the D*tt flight direction instead of the ujtt. 
The polar angle £2 corresponds to the angle £1 but in 
the D*tt rest frame. The cos£i variable is related to 
M 2 (D*tt) whereas the cos £2 is related to M 2 (ujtt). 

Each set of variables (denoted below with the six¬ 
dimensional vector x ) fully defines the kinematics of 
the decay chain, either in the color-favored or the 
color-suppressed channel {36|. The probability den¬ 
sity function (PDF) in the signal region, which is 
the sum of signal and background components, is 
constructed in such a way that the kinematic depen¬ 
dence of the efficiency can be omitted in the mini¬ 
mization [12(1 : 


PDF(x, a) = 


e(x) 


n s T ribkgj 
\M(x,a)\ 2 


.(3) 


£ 


n bkgj 


B j (x) 

e bkg j 


(4) 


where the sum is over the background components 
estimated in the sideband Regions II, III and IV 
(see Fig. [3]), and the efficiencies e s and ebkg j corre¬ 
spond to average signal and background efficiencies, 
respectively, in the signal Region I integrated over 
the phase space. In Eq. 0, a is the vector of pa¬ 
rameters determined from the unbinned likelihood 
fit; n s is the expected number of the signal events 
in the signal Region I distributed according to the 
matrix element squared \M(x, a)| 2 ; e(x) is the recon¬ 
struction efficiency for the B° —> D* + ujtt~ CR events 
in Region I depending on the decay kinematics and 
slowly varying within the scale of resolution of the 
observables; and nbkgj is the expected number of 
background events in the signal Region I distributed 
according to the function Bj(x). We neglect the con¬ 
volution with the resolution function in Eq. 0 due 
to the small invariant mass resolutions (4 MeV/c 2 for 
ujtt and 3 MeV/c 2 for D**) in comparison with the 
resonance widths (more than 150 MeV/c 2 for the p- 
meson-like resonance and more than 25 MeV/c 2 for 
the D** states). 

An unbinned likelihood fit to the B° —> D* + ujtt~ 
phase space is performed to minimize the negative 
log-likelihood function C(a): 

C(a) = - In PDF + 

events 

(?7»s Yjj ^bkg j ^tot) 

2(ntot + o’bkg) 


( 5 ) 
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Figure 6: (color online). Kinematics of a B° — > D* + un decay mediated by an u>7r~ intermediate resonance. The 
diagram in (a) defines two polar angles £1 and /3i and one azimuthal angle tpi. The diagram in (b) defines one polar 
angle Q\ and one azimuthal angle tj> 1 . The direction in (b) corresponds to the vector normal to the u> decay plane. 


where n tot is the total number of events in the sig¬ 
nal Region I and <7bkg is the uncertainty of the total 
number of background events JV ribkgj- The second 
term in Eq. ([5J takes into account our knowledge of 
the background contribution in the signal region. 

The function C(a ) does not incorporate the inter¬ 
ference between the D*Att peaking background and 
the D*ujtt signal. This effect is expected to be small 
(see Section IV.El) . 


A. Background description 

The background components of Eq. a can be 
addressed using the (AE, M (n + n~n° )) scatter plot 
(Fig. a- The combinatorial background with misre- 
constructed w candidates saturates Region IV. The 
B° — > D* + TT + TT~Tr°Tr~ events without uj in the inter¬ 
mediate state can be found in Region III. The com¬ 
binatorial background with a correctly reconstructed 
w falls into Region II. In addition, the SCF events 
lie in all regions. 

We determine the six-dimensional shapes of the 
background PDFs Bj(x ) by performing an unbinned- 
likelihood fit in the sideband regions. For details, see 
Appendix [A] 

The projections on the M 2 (w7t) and M 2 (D*ir) 


variables and the corresponding background fits are 
shown in Fig. [7] The result of the unbinned-like- 
lihood fit in Region IV determining the function 
Biy(x) is shown in Figs. [7] (a) and (b). Figures[7] (c) 
and (d) correspond to Region III. The backgrounds 
in this region are described by the function Bm(x) 
plus a contribution components described by Biy(x). 
In a similar way, Region II includes the background 
components described by Biy(x) and Bu(x); these 
components are shown in Figs. [7] (e) and (f). 

The M 2 (ujit) and M 2 (D*tt) distributions of the 
background in the signal Region I are shown in 
Fig. [8] These distributions are the sum of the SCF 
distribution in Region I obtained from the MC study 
and the distributions describing the backgrounds de¬ 
fined above. The latter distributions are the dif¬ 
ferences between the Bu(x), Bm(x) and Biy(x) 
distributions and the SCF distributions in Regions 
II, III and IV, respectively. Figure [8] illustrates 
the dominant contribution due to the combinato¬ 
rial background with a misreconstructed uj. The 
D* + 7r + 7r _ 7r°7r _ component estimated from Region 
III is also significant in Region I. The combinato¬ 
rial background with a correctly reconstructed uj and 
SCF component obtained from the MC study have 
lower fractions but are also included in the descrip¬ 
tion. 
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M 2 (D* 31) (GeV/c 2 ) 2 



Figure 7: M 2 (ujtv) and M 2 (D*tv) distributions of the B° — > D* + u;7r - candidates in the (AE, M('K + 'K~n 0 )) sideband 
regions (a,b) IV, (c,d) III and (e,f) II. Points with error bars are data; hatched histograms correspond to the 
contribution from dotted histograms represent the component described by the function Bni(a;) in (c) and 

(d) and Bu(x) in (e) and (f); open histograms correspond to the total fit results in Regions III and II. 
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SCF (b) 
co k comb. 

D 4 n 
4 k comb. 





M 2 (co 7i) (GeV/c 2 ) 2 


M 2 (D* ti) (GeV/c 2 ) 2 


Figure 8: (a) M 2 {ujn) and (b) M 2 (D*n) distributions of the background components in the signal region. The 
histograms are stacked on top of each other. 


B. Signal description 

The description of the D* + con~ signal events in 
phase space is based on the study of Ref. {13 ■ Since 
charge conjugation is taken into account, the total 
matrix element is calculated as: 

M = ±±9_ M+ + ^-M_, (6) 

where Q = +1 (—1) for B° ( B° ) decays and M_ dif¬ 
fers from M_|_ by the sign of the P-violating terms. 
Following the isobar model formulation [44j with 
quasi-two-body resonant amplitudes, the matrix el¬ 
ement M± is given by: 

M ± = Y j a R e i * R M R± , ( 7 ) 

R 

where an and <f>n are relative amplitudes and phases 
of the intermediate resonances and R is an index 
numbering all the con and D** resonances. The full 
description of the resonant matrix elements Mr± 
can be found in Appendix [B] The parametrization of 
the form factors used in the matrix elements Ain ± is 
presented in Appendix[Cj The fraction fn of the total 
three-body signal attributed to a particular quasi- 
two-body intermediate state is defined as 

, = / ah\M R± (x)\ 2 p{x)dx 

R J \A4±(x)\ 2 p(x)dx ’ 


where p(x) is the phase space density of events de¬ 
termined from the kinematic conditions of the decay 
{36\. The sum of the fit fractions for all components 
is not necessarily unity because of interference ef¬ 
fects. 

The fraction of resonance R produced in partial 
wave L is determined as 

,l = I \ M R±ix)\ 2 p{x)dx 
J \M R ±(x)\ 2 p(x)dx' 

where ± is the matrix element describing the pro¬ 
duction of resonance R in partial wave L and the sum 
J2 l fk i s unity by definition. 

The observable determined from the amplitude 
analysis is the longitudinal polarization Vn of res¬ 
onance R. This variable is calculated as 

v = _ l g ol 2 _ 

rR \H 0 \ 2 + \H + \ 2 + \H_\ 2 ' ( ’ 

where Ho, H + and F/_ represent three complex helic- 
ity amplitudes which can be expressed via invariant 
and partial wave form factors (see Appendix [Cl). 

C. Fitting the B° —> D* + cj 7r _ signal 

Figure |9] shows the two-dimensional Dalitz distri¬ 
butions in signal Region I and sideband Regions II, 
III and IV. There are 1129 events in the signal region 
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Figure 9: Dalitz distributions of the D* + u>n candidates i 
Region III and (d) sideband Region IV. 

that satisfy all the selection criteria. 

To describe all the features of the Dalitz plot, we 
use the following set of resonances: off-shell p(770)~, 
p(1450)-, £>1(2430)°, £>i(2420)° and £>5(2460)°. A 
CLEO analysis [29( showed the dominance of the 
/?(1450p resonance in this final state. In a BABAR 
study [30|, a D*tt enhancement was observed that 
was interpreted as a £>i(2430)° signal. Our data re¬ 
quire additional resonances. We take into account an 




(a) signal Region I, (b) sideband Region II, (c) sideband 


off-shell p(770)“ contribution, as suggested by the 
e + e _ —» W7r° data [45| . To improve the description, 
we also include the amplitudes of the narrow reso¬ 
nances £>i(2420)° and £>5(2460)°. When both res¬ 
onances are simultaneously included in the matrix 
element rather than just one of them, the statisti¬ 
cal significance of the signal, given by \J2(Cr — Co), 
where Cr (Cq) is the negative log-likelihood value 
with the signal from the resonance R fixed at zero 
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(with the nominal signal yield), increases very sig¬ 
nificantly (> 5cr effect). We also include in the fit 
a SCC contribution with the 6i(1235)“ resonance. 
This contribution has a significance below 3 .Oct and 
we obtain an upper limit for the fraction of the SCC 
in B° -> D*+ ujTT decays. To determine the upper 
limit, we generate pseudoexperiments (see Sec. IV.Dl) . 

The results of the fit are summarized in Table El 
Together with the individual decay fractions /#, we 
show the FCC fraction / p + p /, which represents the 
decay fraction of the coherent sum of the p(770) _ 
and p(1450) _ states. 

We also show the partial wave fractions describing 
the p(1450)“ and D** production in the specific par¬ 
tial waves and the longitudinal polarizations of these 
resonances. We see that the p(1450)“ state is pro¬ 
duced dominantly via S wave, but that D i(2430)° 
production requires approximately equal fractions of 
all partial waves. The partial wave fractions of the 
Di(2420)° and .02(2460)° are not statistically signif¬ 
icant. In our analysis, the longitudinal polarization 
and partial wave fractions of the p(1450)“ are fixed 
in part from the requirement on the relative normal¬ 
izations of the hclicity amplitudes, R±, R -2 and p 2 
(see Appendix O, and can be relaxed in the fit due 
to the free mass and width of the p(1450) - . Large 
longitudinal polarizations of the D** states indicate 
violation of the factorization hypothesis but the sta¬ 
tistical uncertainties are large. 

The final-state interaction phases </>+ and </>_ de¬ 
fined in Appendix [C] are taken into account in the 
description of the p-meson-like states. The fit gives 
a nontrivial value for the </>+ phase. 

One must also consider the statistical errors on the 
fit fractions, partial wave fractions and longitudinal 
polarizations. These errors are determined with a 
pseudoexperiment technique (see Sec. IV.Dl) . 

The masses and widths of all resonances except 
for the p(1450) are fixed at their PDG values |3J. 
Our measurements for the p(1450) shape parameters 
do not contradict previous experimental observations 
[S], although they differ slightly from the CLEO re¬ 
sults [29] . This situation is expected because the 
broad p-meson-like states overlap strongly with each 
other and the Breit-Wigner description is not accu¬ 
rate. 

Mixing between the Z?i(2430)° and Hi (2420)° 
states is expected to be small and is therefore ne¬ 
glected. If we take into account the mixing effect, the 
mixing angles defined in Appendix [C] are found to be 
w = — 0.03±0.02 (stat.) and = — 0.27±0.75 (stat.). 
Within errors, these angles are consistent with the 
previous Belle measurement [12]. 

Figures [TO] and El show the distributions of the 
kinematic variables related to the W7r and D*ir sys¬ 
tems, respectively, for the D* + ujtt~ candidates in the 


signal region. The results of the fit with the nomi¬ 
nal model are superimposed. All plots demonstrate 
a reasonable description of the data by the fit. A 
more detailed comparison is shown in Figs, [l^] and 
m for regions enriched (cos £2 > —0.4) and depleted 
(| cos 61 1 < 0.5) with D** mesons. 

To ensure that our fit results correspond to the 
global minimum, we repeat the signal fit 1000 times 
with randomly selected starting values for the fit pa¬ 
rameters. None of these fits have better likelihoods 
than those presented above. For the nominal fit, 
two local minima are found. One of them, which is 
3.3 ct away from the global minimum, corresponds to 
a very large decay fraction for the p(1450), / p ( 1450 ) = 
(157.3 ±23.1)%, in comparison with the decay frac¬ 
tion for the off-shell p(770), / p ( 770 ) = (87.5± 13.1)%, 
and a relative phase (/> p (i 4 so) = —2.52 ± 0.05. This 
result is inconsistent with the e + e _ data j4f|. For 
the other one, all the fit parameters coincide with the 
values presented in Table nn within their statistical 
errors with the exception of the relative phases in the 
Di (2430)° description: the S- and D -wave phases 
are shifted by 7r/2, whereas the P -wave phase re¬ 
mains unchanged. Since this second local minimum 
is more than 3.5 ct away from the global minimum, 
it is not considered as a second possible solution for 
the final results. 


D. Statistical uncertainties 

The statistical uncertainties on the fractions of 
the different intermediate states as well as the up¬ 
per limit for the b\ (1235) fraction are determined 
with using a frequency method. Another objective 
of this procedure is to estimate how well the nomi¬ 
nal signal model describes the data. Assuming ade¬ 
quate agreement between the data and the nominal 
signal model, we generate 1000 statistically indepen¬ 
dent samples, which are a proper mixture of signal 
and background events distributed according to the 
PDF of Eq. ([4]) . All the major characteristics such as 
the reconstruction efficiency and statistics are taken 
into account. The numbers of signal and background 
events for each pseudoexperiment are generated ac¬ 
cording to distributions by statistics. 

We fit the obtained MC samples and determine the 
fractions of quasi-two-body channels for each sam¬ 
ple. The distributions of these fractions are then fit 
with a Gaussian G(x; p, ct) or bifurcated Gaussian 
(Gaussian with different standard deviation values 
ct 1 and CT 2 on left and right side of the mean value 
p) G(x; p, cti, CT 2 ). The standard deviations, ct or cti 
and <72, are considered as the statistical errors for the 
fractions of the corresponding submode. The 90% 
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Table II: Summary of the fit results to the D* + lutv~ candidates in the signal region. Each column of results corresponds 
to a different signal model. The notations p' = p(1450), D[ = Di(2430), D\ = Di(2420) and D 2 = D 2 ( 2460) are 
used. Quoted uncertainty is statistical only. A£ = C — Co, where C defined in Eq. ([S]) corresponds to the signal model 
for which this variable is calculated and Co is the negative log-likelihood function calculated for the signal model with 
p, p', D'i, D i and D 2 resonances. 


Contribution 

Parameter 

P,p' 

P,p' 

P > p' 

P,p' 

p, pi b i 



D’i 

75i, 75x 

d[,d * 2 

D'uD^DZ 


p(770) - .D* + 

Resonance phase 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 


Resonance fraction, % 

65.7 ±10.2 

65.2 ± 11.7 

61.6 ±11.9 

64.2 ± 10.7 

60.6 ± 12.1 

p(1450) _ D* + 

Resonance phase 

2.63 ±0.11 

2.55 ±0.11 

2.62 ±0.11 

2.56 ±0.12 

2.54 ±0.11 


Resonance coupling 

0±8t°;° 2 

n i q+ 0-02 

U--L°-0.05 

0 20 -1 " 0 ' 03 
U.ZU—o.06 

n i q+ 0-02 
U.lo_ 0 06 

n 10+0.03 
U.ly_ 0 .06 


Mass, MeV/c 2 

1549 ± 22 

1546 ± 23 

1543 ± 23 

1544 ± 22 

1540 ± 22 


Width, MeV/c 2 

3031“ 

3051“ 

316±“ 

3031“ 

3021“ 


Ri 

1.40 (fixed) 

1.40 (fixed) 

1.40 (fixed) 

1.40 (fixed) 

1.40 (fixed) 


R2 

0.87 (fixed) 

0.87 (fixed) 

0.87 (fixed) 

0.87 (fixed) 

0.87 (fixed) 


p 2 

Resonance fraction, % 
S- wave fraction, % 

0.79 (fixed) 
46.71H9 
76.911 2 

0.79 (fixed) 
44.5l® 2 9 0 
75.311? 

0.79 (fixed) 
50.41“;? 
76.71?;? 

0.79 (fixed) 
46.3l? 3 ° 4 
75.11?;? 

0.79 (fixed) 
47.5l? 2 3 1 
75.3±f;g 


P- wave fraction, % 

12.0 ±0.7 

12.8 ±1.1 

12.2 ±0.9 

12.9 ±0.9 

12.7 ±0.8 


D- wave fraction, % 

11.0 ±0.4 

11.8 ±0.6 

11.0 ±0.6 

11.9 ±0.5 

11.8 ±0.6 


4>+ phase 

0.66 ±0.33 

0.86 ± 0.30 

0.67 ±0.37 

0.87 ±0.29 

0.85 ± 0.31 


<f>- phase 

-0.14 ±0.17 

-0.02 ±0.15 

-0.15 ±0.19 

-0.02 ±0.13 

-0.02 ±0.15 


Long, polarization, % 

66.4 ±0.6 

66.5 ± 0.6 

66.5 ±0.6 

66.5 ± 0.6 

66.6 ±0.6 


FCC fraction, % 

79.1 ±2.5 

82.6 ±2.4 

79.0 ±2.4 

82.2 ±2.2 

81.6 ±2.3 

Di(2430) u a; 

Resonance phase 

0.91 ±0.26 

1.03 ±0.28 

1.11 ±0.29 

1.24 ±0.28 

1.27 ±0.35 


S- wave phase 

0.26 ±0.20 

0.19 ±0.23 

0.14 ±0.23 

-0.05 ±0.25 

-0.09 ± 0.26 


P- wave phase 

2.71 ±0.21 

2.41 ±0.27 

2.56 ± 0.24 

2.24 ±0.29 

2.23 ±0.32 


Resonance fraction, % 

13.6 ±2.1 

11.2 ± 1.8 

12.6 ±1.8 

10.8 ±1.8 

11.6 ±2.0 


S-wave fraction, % 

29.7 ±8.6 

33.6 ± 9.5 

35.8 ± 10.1 

38.9 ± 10.8 

38.9 ± 10.5 


P- wave fraction, % 

37.0 ±8.6 

34.1 ±9.2 

34.0 ± 8.9 

33.1 ±9.5 

29.1 ±9.1 


D- wave fraction, % 

33.5 ±8.8 

32.6 ±9.2 

30.5 ± 9.2 

28.3 ±8.9 

32.2 ±9.2 


Long, polarization, % 

60.9 ±8.2 

63.4 ±8.9 

63.0 ±8.2 

63.0 ±9.1 

67.6 ±9.2 

Di(2420) u a; 

Resonance phase 


1.92 ±0.34 


2.12 ±0.34 

2.16 ±0.42 


S-wave phase 


-0.06 ±0.34 


-0.07 ±0.43 

-0.10 ±0.43 


P- wave phase 


0.04 ±0.41 


-0.25 ± 0.46 

-0.24 ±0.49 


Resonance fraction, % 


3.7 ± 1.1 


2.9 ±0.8 

2.8 ±0.8 


S'-wave fraction, % 


35.6 ± 13.2 


34.0 ± 13.4 

35.8 ± 13.0 


P -wave fraction, % 


36.6 ± 11.8 


31.2 ± 11.4 

30.3 ±11.0 


75-wave fraction, % 


27.9 ± 11.0 


34.9 ± 13.4 

34.0 ± 13.1 


Long, polarization, % 


60.2 ± 12.0 


67.1 ± 11.7 

67.4 ± 16.1 

D 2 (2460)°u; 

Resonance phase 



1.69 ±0.57 

2.31 ±0.50 

2.39 ± 0.42 


P-wave phase 



-0.67 ±0.54 

-0.77 ±0.62 

-0.84 ±0.52 


75- wave phase 



-1.10 ±0.71 

-1.85 ±0.59 

-1.96 ±0.58 


Resonance fraction, % 



2.1 ±0.7 

1.8 ±0.6 

1.8 ±0.6 


P- wave fraction, % 



34.3 ± 16.6 

29.5 ±16.9 

30.0 ± 16.7 


75-wave fraction, % 



45.7 ± 17.4 

40.2 ± 17.7 

38.2 ±17.3 


P-wave fraction, % 



19.4 ± 15.8 

29.4 ± 19.3 

31.1 ± 19.2 


Long, polarization, % 



74.1 ± 16.5 

76.0l? 8 5 3 

74.7 ± 16.1 

6i(1235) _ L>*+ 

Resonance phase 





0.52 ±0.42 


Resonance fraction, % 





< 3.1 (90% C.L.) 


AC 

±33.3 

±12.9 

±16.4 

0 

-2.4 


Variation, <r 

8.2 

5.1 

5.7 

0 

2.2 
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Figure 10: (color online). Distribution of the six unr variables for D*cun candidates in the signal region (points with 
error bars). The histograms represent the results of the fit (black), including the following components: p(770) (cyan), 
p(1450) (red), p( 770) and p(1450) together (red dashed), Di(2430)° (green), Di(2420)° (blue), DJ(2460)° (magenta) 
and background (hatched). 
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Figure 11: (color online). Distribution of the six D* tt variables for D*u>n candidates in the signal region (points with 
error bars). The histograms represent the results of the fit (black), including the following components: p(770) (cyan), 
p(1450) (red), p( 770) and p(1450) together (red dashed), Di(2430)° (green), Di(2420)° (blue), DJ(2460)° (magenta) 
and background (hatched). 
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Figure 12: (color online). Distribution of three 0J7T variables for Z)* + a;7r _ candidates in two different subregions of 
the signal region (points with error bars), defined by cos £2 > —0.4 or | cos(9i| > 0.5 ( D ** enriched) and | cos(9i| < 0.5 
(D** depleted). The histograms represent the results of the fit (black), including the following components: p(770) 
(cyan), p(1450) (red), p( 770) and p(1450) together (red dashed), Di(2430)° (green), Di(2420)° (blue), D5(2460)° 
(magenta) and background (hatched). 
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Figure 13: (color online). Distribution of three D*ix variables for D* + lo-k~ candidates in two different subregions of 
the signal region (points with error bars), defined by cos £2 > —0.4 ( D ** enriched) and |cos#i| < 0.5 ( D ** depleted). 
The histograms represent the results of the fit (black), including the following components: p(770) (cyan), p(1450) 
(red), p(770) and p(1450) together (red dashed), Di(2430)° (green), Di(2420)° (blue), DJ(2460)° (magenta) and 
background (hatched). 
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confidence level upper limit for the 6i(1235) _ con¬ 
tribution is obtained directly from the distribution 
of the &i(1235)“ decay fractions in the pseudoexper¬ 
iments. 

To measure the goodness of the fit, we utilize two 
different approaches. The first operates with a mixed 
sample [461] combining the experimental data sample 
and pseudoexperiments with ten times higher statis¬ 
tics than in the experiment. This method allows 
one to estimate the consistency of the nominal signal 
model and data in the multidimensional amplitude 
analysis with the small data sample when the \ 2 
method with binning in the multidimensional phase 
space is not valid._Following the algorithm described 
in detail in Ref. [Tti] . we conclude that our nominal 
model and the data are consistent at 49% confidence 
level. 

For the second technique, we define two y 2 vari¬ 
ables calculated in the ujtt and D*n bases, respec¬ 
tively. For each pair of kinematic variables j and k 
describing the W7r (or D*n) production, we consider 
10 x 10 = 100 two-dimensional bins and compute 


2 _ Y' — Nobs i) 2 

N\~- 

-i iY obst 

i=l 


(ii) 


In Eq. ED, Nati is the expected number of events 
in bin i based on the PDF of Eq. 0 and N ohsi is 
the number of observed events in that bin. Then 
we obtain the total \ 2 as the sum of y 2 k over all 
possible pairs of variables j and k. In 90% (78%) of 
the pseudoexperiments this % 2 , calculated with the 
ujTT (D*tt) variables, has a value smaller than in the 
data, indicating an acceptable fit quality. 


E. Systematic uncertainties 

Two types of uncertainties are considered besides 
the statistical errors. These are systematic and 
model uncertainties. 

The systematic uncertainty comes from the back¬ 
ground description and the efficiency of the selection 
requirements. To estimate the uncertainty in the 
parametrization of the distribution of background 
events, we use two alternative parametrizations. The 
first is determined in terms of the D** production 
variables instead of the ujit basis used in the nomi¬ 
nal fit. In the second, we use alternative background 
functions: a sum of Legendre polynomials instead 
of a sum of exponential functions used in the nom¬ 
inal fit and alternative correlation functions. The 
full parametrization for the nominal background fit 
is presented in Appendix [A] The uncertainty related 
to the efficiency of the definition of the selection re¬ 


quirements is dominated by the variation of the sig¬ 
nal region in the plane. To es¬ 

timate this uncertainty, we modify the signal region 
shape from the rectangle to an ellipse, taking into ac¬ 
count the correlation between A E and 
This modification increases the signal-to-background 
ratio by a factor of about 1.5. The contributions to 
the uncertainty from the background description and 
the reconstruction efficiency are added in quadrature 
to obtain the overall systematic uncertainty. 

The uncertainties in the parametrization of the 
signal matrix element determine the model error. 
There are three sets of such uncertainties. The first is 
related to the number of contributions to the matrix 
element. We include an additional nonsignificant 
6i(1235)~ amplitude described in Appendix IBl and 
then try several fits: first, we modify the 6i(1235)“ 
rnodel by removing the D -wave contribution to the 
decay &i(1235)“ —> then we consider the rel¬ 

ative helicity phases (f>± as free parameters during 
the fit, independent of the values of the helicity 
phases defined in the p( 770) and p(1450) ampli¬ 
tudes. The model error due to the &i(1235) con¬ 
tribution is assigned as the maximum difference be¬ 
tween the values obtained from these fits and the 
nominal one. Furthermore, we include in the sig¬ 
nal model the contributions from p(1700) _ , off-shell 
D° resonances and S -wave nonresonant amplitudes. 
All of them are nonsignificant. The second set of 
errors arises due to the assumption of the signal 
shape. We take into account the mixing effect be¬ 
tween the Z?i(2430)° and Z?i(2420)° states. More¬ 
over, we modify the transition form factors in the 
matrix element: we substitute the effective form fac¬ 
tor A(q 2 ), describing the p(770) _ —> ojn~ transition 
for the P -wave Blatt-Weisskopf factor Bp(q 2 ) (see 
Appendix [Bj and we modify the shape of the Isgur- 
Wise function h{w ) describing the production of the 
p- meson-like states (see Appendix [C]) . For the lat¬ 
ter, we apply the parametrization that corresponds 
to the requirements of analyticity and is used in the 
BABAR B° —> D* + e~v e analysis [47j| . The third set 
of errors is related to the model parameters that are 
fixed in the fit. We vary the mass and the width 
of each resonance [except the p(1450)] within their 
known PDG uncertainties Q. We also vary the pa¬ 
rameters R± , i ?2 and p 2 of the invariant form fac¬ 
tors describing the p-meson-like amplitudes (see Ap¬ 
pendix [C]) within their uncertainties obtained by the 
BABAR collaboration J47|. Moreover, we vary the 
parameter r = 1.6 (GeV/c) -1 used in the Blatt- 
Weisskopf factors and the form factor A(q 2 ) (see Ap- 
pendix IBl) in the range from 0.8 to 2.5 (GeV/c) -1 . 

The total model error is obtained by adding all 
model errors in quadrature. The sources of system¬ 
atic and dominant model uncertainties that affect 
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Table III: Summary of systematic and dominant model uncertainties in the parameters of amplitude analysis. 


Contribution 

Parameter 

Systematic 



Dominant model 




uncertainties 


uncertainties 




Background 

Signal 


r 

Mixing 



description 

region 

form factor 

parameter 

effect 

p ( 770 )" D * + 

Resonance fraction, % 

+ 0.8 

- 0.1 

± 4.9 

± 6.1 

+F 77 

- 24.0 

- 0.9 

p ( 1450 )- L >*+ 

Resonance phase 

+ 0.01 

- 0.02 

± 0.07 

± 0.24 

+ 0.22 

- 0.17 

0.00 


Resonance coupling 

± 0.01 

- 0.02 

± 0.08 

+ 0.10 

- 0.01 

± 0.01 


Mass, MeV/c 2 

±1 

±11 

-17 

+ 1 
-42 

0 


Width, MeV/c 2 

+2 

-4 

±3 

±69 

+55 

-6 

±2 


Resonance fraction, % 

± 1.9 

- 4.4 

± 9.9 

+ 17.4 
- 0.8 

± 0.7 


</>+ phase 

± 0.05 

± 0.07 

± 0.06 

± 0.06 

0.00 


</>_ phase 

- 0.02 

± 0.02 

± 0.05 

± 0.03 

- 0.01 


FCC fraction, % 

- 0.2 

- 3.6 

± 0.3 

+ 0.3 

- 1.8 

- 0.5 

Di ( 2430 ) u cu 

Resonance phase 

- 0.07 

± 0.18 

- 0.29 

+ 0.39 

- 0.32 

± 0.03 


S-wave phase 

± 0.04 

- 0.26 

± 0.04 

+ 0.04 

- 0.02 

- 0.05 


P -wave phase 

+ 0.03 

- 0.04 

- 0.26 

± 0.13 

+ 0.08 

- 0.05 

- 0.04 


Resonance fraction, % 

± 0.1 

± 2.7 

- 0.4 

+ 1-1 
- 0.2 

± 1.3 


S'-wave fraction, % 

+ 4.2 

- 0.7 

± 0.9 

- 0.3 

- 1.0 

± 1.2 


P -wave fraction, % 

+ 1.2 
- 5.5 

± 2.1 

- 0.3 

+ 2.9 

- 0.1 

± 0.8 


D -wave fraction, % 

- 0.8 

± 3.0 

± 0.5 

+ 0.2 

- 2.0 

- 2.1 


Long, polarization, % 

+ 4.6 

- 1.2 

- 4.4 

± 0.4 

+ 0.6 

- 3.5 

- 1.8 

Di ( 2420 ) u u ; 

Resonance phase 

+ 0.08 

- 0.03 

± 0.08 

- 0.23 

+ 0.32 

- 0.27 

± 0.05 


S- wave phase 

+ 0.03 

- 0.17 

± 0.09 

± 0.11 

+ 0.05 

- 0.07 

± 0.04 


P -wave phase 

± 0.07 

- 0.37 

± 0.02 

+ 0.02 

- 0.04 

± 0.03 


Resonance fraction, % 

± 0.2 

± 0.4 

- 0.2 

+ 0.0 

- 0.1 

± 0.5 


Long, polarization, % 

- 3.7 

- 2.0 

- 0.9 

± 2.0 

- 2.8 

D 2 ( 2460 ) u tu 

Resonance phase 

± 0.03 

- 0.12 

- 0.24 

± 0.30 

± 0.03 


P- wave phase 

+ 0.02 

- 0.11 

- 0.10 

± 0.04 

+ 0.02 

- 0.12 

- 0.04 


D- wave phase 

+ 0.01 

- 0.06 

- 0.37 

± 0.08 

± 0.07 

- 0.08 


Resonance fraction, % 

+ 0.0 

- 0.1 

0.0 

0.0 

+ 0.0 

- 0.1 

± 0.1 


Long, polarization, % 

+ 0.2 

- 2.0 

± 2.0 

± 1.5 

+ 1.4 
- 0.3 

- 1.5 


the results of the amplitude analysis are summarized 
in Table Hill 

To account for the systematic and model uncer¬ 
tainties in the upper limit of the &i(1235) _ , we de¬ 
termine the &i(1235) - contribution with all above 
described sources of errors [including the 6i(1235) 
mass and width variation] and use the largest value 
to evaluate the upper limit. The main effect is due to 
the removal of the D wave in the fri(1235)“ -± ojtt~ 
decay. 

An additional effect appears due to the interfer¬ 
ence between Z?* + 7r + 7r”7r°7r” background events and 
D* + um~ signal events. Figure [8] (a) shows that 
most of the D*Air events lie in the range M(4n) < 
2GeV/c 2 . The investigation of e + e~ annihilation 
into a 4-7r system [48j] at these energies as well as 
the study of the resonant structure in the decay 
r —> 37r7r°jz T [l!| demonstrate the dominance of 
the ai(1260)7r and ojtt intermediate states. We as¬ 
sume that our D*4ir background is also dominated 


by ai(1260)7r production. In such a case, the inter¬ 
ference with the ujtt system should be negligible. 


VI. DISCUSSION AND CONCLUSION 

This analysis is devoted to the study of the three- 
body B° —> D* + ojn^ decay. We obtain the total 
branching fraction 

B = (2.31 ± 0.11 (stat.) ± 0.14 (syst.)) x lO" 3 , 

consistent within errors with the CLEO [29] and 
BABAR [30] measurements but with a slightly 
smaller central value. 

A full amplitude analysis of the final state has been 
performed. A summary of the results with system¬ 
atic and model uncertainties on parameters and sta¬ 
tistical significances of resonant contributions is pre¬ 
sented in Table HVl This is the first consistent study 
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Table IV: Summary of the final results of the B° —> D* + con~ amplitude analysis. The first error is statistical, the 
second is systematic and the third is the model error. The statistical significance, taking into account systematic 
effects, is given by \J 2(A£), where AC is the difference between the negative log-likelihood values with the signal 
from resonance fixed at zero and the nominal signal yield. 


Contribution 

Parameter 

Value 

Significance 


Total branching fraction, ICC 3 

2.31 ±0.11 ±0.14 



FCC branching fraction, 10 -3 

i QQ r» i i + 0.11 + 0.02 

i.yu m u.ii_ 0 . 13 _o.o 6 

29.8cr 


SCC branching fraction, 10 ~ 4 

< 0.7 (90% C.L.) 


p(770)".D* + 

Resonance phase 

0 (fixed) 



Resonance coupling 

1 (fixed) 



Resonance branching fraction, 10 -3 

1.48 ± 0.27i°;J® j :°; 24 

10.5(7 

p(1450)”D*+ 

Resonance phase 

Resonance coupling 

Mass, MeV/c 2 

Width, MeV/c 2 

2.56 ±0.12^;™!) 

r\ -| O+ 0 . 02 -|- 0 . 00 + 0.10 

U.lO_ 0 .06-0.02-0.01 
1544±22± 31 t 4 6 

orjo+31 +3 +69 

OUO_52_4_6 



Resonance branching fraction, 10 -3 

n7 +0.15+0.06+0.40 
l-U/_o.3i_o.13-0.02 

15.0(7 


4>+ phase 

0.87 ± 0.29±g'o7 ± 0.06 



<f>- phase 

-0.02 ±0.13 ±0.02 ±0.05 


Di(2430)°(u 

Resonance phase 

S- wave phase 

P- wave phase 

1.24 ± O. 28 +Q Q 7 +Q 32 
-0.05±0.25t°;™7 

2.24 ± 0.29t°;°2ito.ol 



Resonance branching fraction, 10 -4 

2.5±0.4t°™ 

8 . 6 cr 


S-wave fraction, % 

38.9 ± 10.8lo 3 !] ? 



P- wave fraction, % 

33.1±9.5;+ 4 t 3 ;° 



D- wave fraction, % 

28.3 ±8.9 1 3 ;°±!;® 



Long, polarization, % 

63.0 ± 9.1 ±4.6^3^ 


Di(2420)°o; 

Resonance phase 

S- wave phase 

2 12 + 0 Q 4 + u ll + u - 3a 

-E u ‘ O 4 : - 0 . 03 - 0.27 

_n Q 7 4 - n 40 + 0 . 09 + 0.12 

U.UI zl U.^O—o.l7—0.08 



P- wave phase 

-0.25 ± 0.46tg;^ ± 0.04 



Resonance branching fraction, 10 -4 

0.7 ± 0.2l°;J ± 0.1 

5.5cr 


Long, polarization, % 

67.1±11.7l^j: 2 j 


DJ(2460)°w 

Resonance phase 

P- wave phase 

D- wave phase 

2.31 ±0.50+^ ±0.11 

-0.77 ±0.62+™]™ 

-1.85 ± 0.59t°° 4 t°°? 



Resonance branching fraction, 10 -4 

0.4 ± 0.1±g ? ±0.1 

5.0cr 


Long, polarization, % 

76.0± 48 5 3 ± 2.0±1;® 



of the p(770) and p(1450) states in R-meson decays. 
Large signals correspond to the off-shcll p(770)~ me¬ 
son and p(1450) - resonance with significances of 
10.5cr and 15. Oct calculated from the negative log- 
likelihood values and taking into account systematic 
effects. However, model uncertainties are of about 
40%. There is no accurate description yet of the 
shape of the p(1450) resonance. This leads to an 
ambiguity in discriminating the p-meson-like states 
and to large model errors in the definition of their 
resonance branching fractions. Nevertheless, the co¬ 
herent contribution of these resonances is determined 
with smaller model uncertainties. The statistical sig¬ 
nificance of this fraction is 29.8(7. This combined 
decay fraction gives a dominant contribution to the 
total branching fraction. 


We also measure the relative coupling and the rela¬ 
tive phase between the p-meson-like states. Neglect¬ 
ing final state interactions for p-meson-like produc¬ 
tion, we can compare this production with the e + e^ 
SND data j5f| • In the SND analysis j45] as well as in 
our analysis, a small change of the resonance shape 
leads to significant shifts in the fitted resonance pa¬ 
rameters. However, within the isotopic invariance 
and CVC fitted p-meson-like resonance parameters 
are compatible to those observed in processes pro¬ 
ceeding through a virtual photon in e + e - collisions 

III, Ei. 

The phase difference between the p(770)~ and 
p(1450)“ amplitudes is observed to be close to tt as 
predicted in Ref. [51]. In the frame of our signal 
model, we measure the p(1450) _ mass and width. 
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Our measurements also show evidence for nontriv¬ 
ial final-state interaction phases in the helicity am¬ 
plitudes of the p-meson-like states, off-shell p(770)^ 
and p(1450)“, with a significance of 3.3er. Such ef¬ 
fect is observed within the validity of the factoriza¬ 
tion at relatively low q 2 . We restrict the descrip¬ 
tion of these resonances by the requirement that 
the helicity phases cf> + and <f>_ in the amplitude of 
the off-slrell p(770)“ are equal to the corresponding 
phases in the amplitude of the p(1450) _ . Similar 
phases were measured by the CLEO collaboration 
in B —> D* p — > decays Our results as 

well as those of CLEO show that > 4>~, although 
statistical uncertainties are large. 

In addition to the p-meson-like states, the 
&i(1235)“ resonance could be produced as a pos¬ 
sible intermediate state in the color-favored chan¬ 
nel. Such a contribution is generated by SCC and 
is expected to vanish in the limit of perfect isospin 
symmetry. Our measurements do not require any 
SCC contribution and an upper limit for the product 
of branching fractions of B(B° —> £>* + &i(1235) _ ) x 
£>(&i(1235)“ —» W7T - ) has been obtained. This result 
is the first search for SCC in £?-meson decays. 

Color-suppressed decays B° — > £>i(2430)°u; and 
B° —> £>i(2420)°w are observed in our study with 
significances of 8.6er and 5.5tr, taking into account 
systematic effects. The measurements show the rel¬ 
ative dominance of the broad Di (2430)° production 
in comparison with the narrow £>i(2420)°. Heavy 
quark symmetry predicts the absence of Di (2420)° 
signal in the limit A qcd /me ~> 0 [lf|, where A qcd 
is the QCD scale and m c is the mass of the c quark. 
The production of the £>i(2420)° state can be ex¬ 
plained by finite corrections of order Aqc_d/to c in 
the £>i(2420)° production. Moreover, the dominance 
of broad resonances in the color-suppressed channel 
can result in comparable production of the broad and 
narrow states in the B~ —> D**°n~ decays [l3 |. 

The nonfactorizable B° —> £>2(2460)°w decay 
has been observed with a statistical significance of 
5.0cr. In SCET theory 0, the equality of branch¬ 
ing fractions and strong phases in the decays B° —» 
D*(2460)°M and B° £>i(2420)°M, where M = 
it, p, K or M = K* with longitudinal polarization, is 
predicted. Our result with M = w also does not con¬ 
tradict this prediction. However, our errors on the 
resonance branching fractions and phases are large. 

In our analysis, we obtain the partial wave frac¬ 
tions for the intermediate resonances. For the p- 
meson-like states, we fix the relative normalizations 
R\ 1 i ?2 and parameter p 2 in the helicity amplitudes 
at values obtained from the semileptonic B —»• D*lv 
analysis |47|. These normalizations determine the 
relative partial wave fractions with the dominant 
S- wave production (see Table [II]). Another effect 


takes place for the ££(2430)° resonance: all partial 
waves— S, P and D waves—have close probabilities 
of about 30%. A similar tendency is observed in the 
£>i(2420)° and £>2(2460)° production. However, the 
statistical accuracy is not sufficient to obtain signif¬ 
icant numerical values (see Table HT1). 

We also measure for the first time the longitudinal 
polarization of the w in case of D** production. The 
results have large errors, but they imply nontrivial 
nonfactorizable QCD effects in the color-suppressed 
channel [53j] and can be compared with the measure¬ 
ment of the polarization in the decay B° —> D*°u> 
[26| . Vd * = (66.5 ± 5.0)%. All these polarization 
results, except for the £>2(2460)°, show significant 
deviations from unity. The £>2(2460)° result should 
be considered separately because this tensor state is 
generated only due to nonfactorizable contributions. 
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APPENDIX A: BACKGROUND 
DESCRIPTION 


In this section, we present the procedure of the 
background description in the (AE, M(Tr + Tr~n 0 )) 
plane. 

To describe the combinatorial background with¬ 
out real ui events, we use Region IV, which includes 
a SCF component. In this region, a negative log- 
likelihood function to be minimized is given by 


Av(«) 


~H ln 


/ Riv(q,£)e(f) \ 

VScR Av(<?, x) ) 


(A.l) 


where the sum E(rv) over the events in Region 
IV, the sum Ecr calculated over B° D* + ujtt~ 
CR events, which are uniformly generated over the 
phase space and then reconstructed in Region I with 
the above-described selection procedure, and e(x) is 
the reconstruction efficiency for the B° —> D* + ojtt~ 
CR events in Region I. 

After the estimation of Av, the minimization pro¬ 
cedure is performed for events in Region III. In addi¬ 
tion to the events described by the function Av> this 
region includes B° —> D* + TT + Tr~TT°n~ events with¬ 
out oj in the intermediate state and another SCF 
component. The negative log-likelihood function in 
Region III is 


An (b) 


5 > 


/ An-Av Biy(a, x)e(x) 

V Siv A) M Ecr Av(«> %) 

_I 


A AiiAv \ Bm(b,x)e(x) \ , 

l ^vaWEcr^iiiM;’ ' 


where Aii (Siv) is the size of Region III (IV), JVm 
(Niv) is the number of events in Region III (IV), the 
sum E(in) is over the events in Region III and the 
sum Ecr calculated over CR events. The vector a 
in the function Byv is obtained from Region IV and 
fixed in Region III. The vector b is free in Region 

III. 

A similar procedure is performed in Region II. In 
addition to the events described by the function Rrv, 
Region II includes the combinatorial background 
with a real w and another SCF component. The 
shape function Bn describes these events together 
with the additional SCF as in to Region III. The 
minimization function Ai(c) is similar to An(6). 

Functions Bn, An and Ay describe specific back¬ 
ground components defined above and SCF events in 
Regions II, III and IV. All these background contri¬ 
butions are present in signal Region I. However, the 
signal region also includes additional SCF in compar¬ 
ison with the SCF level obtained from the sideband 
regions. This additional SCF component is deter¬ 
mined in MC simulation that shows the same phase 
space distribution in x for all SCF events in each re¬ 
gion of the (A E, M(n + Tr~n 0 )) plane. We can repeat 
a fit in any of the sideband Regions j = II, III or 

IV, taking into account this contribution, and thus 
obtain more precisely the function Bj, which now 
describes this SCF and is used further in the signal 
fit. We choose Region II with the function Ai- 

We use the following empirical parametrization to 
describe the distribution of background events: 


Bj(x) = Fi(M 2 (um), cos A) A (cos A) x 

A(A)A(cos /3i)A(A), (A. 3) 


where the function Fi(A/ 2 (u;7r), cos£i) describes the 
correlation between the M 2 (un) and cos£i variables: 


A (M 2 (w7r), cos £ i) = (gCiv^cosU + _ cos £i) 3 e C4A + c 5 y/A( 1 + cos£i)e C6A ), (A.4) 


and the functions F 2 (cos0i), A^i), A (cos A)> A(V’i) describe one-dinrensional projections of the other 
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variables: 

F 2 ( cosfli) = e C7COs()l +c 8 e C9COsf>1 , 

^ 3 (^ 1 ) =1 + ci 0 sin 2 (^i), 

F 4 ( cos ft) = e CllCOS/31 +ci 2 e Cl3COS/31 , 

^ 5 (^ 1 ) =1 +ci4sin 2 (^i). (A.5) 

Here, d are free parameters, A = — Mq (ujtt) 

and the lower boundary Mq(ujtt) = 0.7 (GeV/c 2 ) 2 
differs from the kinematic limit (m w -\-mA ) 2 because 
the ui invariant mass is not constrained to its nominal 
value. 


APPENDIX B: RESONANT AMPLITUDES 


In this section, we present all resonant amplitudes 
used in the fit. The notations p 2 = M 2 (tt + n~ n°), 
where 7r + 7r _ 7r° is the ui decay product system, and 
q 2 = M 2 (unA) ( q 2 = Ad 2 {D* ± 71 ^)) for the p-meson- 
like (D**) production in the B°(B°) —» D* ± Ld'ir zf: de¬ 
cay, are used. The magnitudes of the three-momenta 
of the oj decay product system and D * ± in the uiit t 
and D^iA rest frames are denoted as p 3 ^ and 
Pd*, respectively. The magnitude of the w three- 
momentum in the coir ^ rest frame when M ( 7 r + 7 r _ 7 r°) 
is equal to the w nominal mass is denoted as p u . 
The magnitude of the oj (D*±) three-momentum 
in the iotA (D* ± tt^ 1 ) rest frame, when M(oj 

is equal to the nominal mass of the p- 
meson-like (D**) resonance and M( 7 r + 7 r - 7 r°) is equal 
to the ui nominal mass, is denoted as p o, w (jPn. p*) - 
The Blatt-Weisskopf penetration factors Bp(p) 541 
used in the resonant matrix element description are 
defined for L = S, P, D and F partial waves as 


Bs(p) 

Bp ip) 

Bd{p) 

Bf(p) 


l + x l 

1 + x 2 ’ 


V 2 -3) 2 +9x 2 
(. x 2 — 3) 2 + 9a: 2 ’ 


I Xq(xq — 15) 2 + 9(2xq — 5) 2 


x 2 {x 2 - 15) 2 + 9(2a: 2 - 5) 2 


,(B.l) 


where x = rp, xo = rpo, r = 1.6(GeV/c) _1 is the 
hadron radius and p and po are the magnitudes of 
the daughter particle three-momenta in the mother 
particle rest frame for the case when the resonance 
invariant mass squared is equal to q 2 and the nominal 
mass squared, respectively. 

B°(B°) -a D* ± p(770) =F -A D^wttF 


Since the off-shell p(770) T has J p = 1 _ , the 
pair D * ± and ojtt ^ can be produced in three par¬ 
tial waves: 5, P and D. S and D waves violate 
C- and P- parities and have the additional phase 
7 t/ 2 in comparison with a P wave. The u> and tA 
pair is produced in a P wave via the strong decay 
p(770) T -A <jJlA. 

The resonance matrix element M p ± describing the 
p(770)^ contribution in the B°(B°) -A D^uin^ de¬ 
cay is 


M p ± = 


V¥p3nA(p 37I ) 

DM 2 ) 


fp{q 2 )App ± 


ifs(q 2 )Asp ± ifD(q 2 )A.DP 


(B.2) 


where A(p 3v ) is the effective form factor describing 
the p T -A urn* transition, fs(q 2 ), fp(q 2 ) and /d(<? 2 ) 
are the partial wave form factors obtained in Ap¬ 
pendix [C] Asp, App and Adp are the angular de¬ 
pendencies shown in Table IVl that correspond to the 
definite partial waves in the B°(B°) -A D* ± p(770) p 
and p{ 770) =F -A un A decays, and D p (q 2 ) is the Breit- 
Wigner (BW) denominator, describing the p(770)^ 
shape: 

DM 2 ) = q 2 -m 2 + im p TM 2 )- (B.3) 

Here, m p is the p( 770) T mass and r p (g 2 ) is the q 2 - 
dependent width. 

The form factor A(pM restricts an overly rapid 
growth of the matrix element of the decay (A -A unA 
with p 3 ^ and is chosen as (H) 

= TTE^F' (B - 4) 

The width T p (g 2 ) for events with q 2 > {in u + 
TOtt) 2 , where m^ is the mass of the charged pion, 
is parametrized as [55| 


r P (<? 2 ) = 


p k * b 2 p {k)t p + 9 “ p * A2{jPui) ^ - 3 


\[<f ^ 0 


12tt 


m n 


-P^ 


(B.5) 

Here, g upn is a coupling constant, which is equal to 
16 (GeV/c 2 ) -1 (56), k w is the magnitude of the mo¬ 
mentum of the tA in the p^ -A n^n 0 decay com¬ 
puted in the p{ 770) T rest frame, is the same 
magnitude, when i/g 2 = m p = (775 ± l)MeV/c 2 , 
and T p = (149 ± 1) MeV/c 2 is the p(770) T width Q. 
The first term in Eq. (IB.51) corresponds to the dom¬ 
inant p(770) T decay mode to the r° system and 
the second term describes the wp7r interaction. For 
events with q 2 < (m u + m,) 2 , we use T p (q 2 ) = T p . 

The magnitude and phase, corresponding to this 
resonant amplitude, are fixed at values 1 and 0, re- 
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Table V: Angular dependencies corresponding to the u>7r p 
quantum numbers J p = 1~. Li (L 2 ) is the relative or¬ 
bital angular momentum between the D * ± and (w 
and 7r"^). The notations c a = cos a and s a = sin a are 
used. The angles 9, <j>, /?, ip, £ correspond to the 
angular basis. 


Li 

l 2 

1/2 

S 

P 

—seScjjCps^ + sec^sps^, — ses^spc^c^ 

P 

P 

SQScfrSpS^Ct: H - SgC(j)SpCijj 

D 

P 

2sos$cpS£ + sec^spsp — ses^spcs/jC^ 


spectively. The free parameters are the relative he- 
licity phases <p± in the form factors fs(q 2 ), fp{q 2 ) 
and f D (q 2 )- 

B°(B°) ->• D* ± p(1450) :F -t D* ± w7r =F 

The resonant matrix element corresponding to the 
p(1450) T intermediate state has a form similar to 
Eq. (IB.2[) except for the form factor A(p 3 /) and the 
width T p (q 2 ). Since the ^(MSO)^ is on-shell, we use 
the Blatt-Weisskopf form factor Bp(p 37I ) instead of 
A(p 3 n) ..55||. The width r p ( 145 o) (q 2 ) for events with 
q 2 > (m w + m,r) 2 is parametrized as [55| 


The free parameters for the amplitude 

obtained from the fit are the relative magnitude and 
phase, the mass and width of the ^(ldSO)^, and the 
helicity phases <p±, which are the same as in the 
p(770) p amplitude. 

B°(B°) ->• D* ± bi(1235) =F 

The 6i(1235) =F resonance has quantum numbers 
J p = 1 + . As such, its wave function has an addi¬ 
tional phase 7 t/ 2. The resonant matrix element is 
written as 


M bl ± = 


mb-tBsipzir) ± f P {q )Aps + 


D bl (q 2 ) 

ifs{q 2 )Ass + ifD(q 2 )ADs^j - 
ODse^ DS Pl(p 37 r)B D (p 3 n) X 
± fp{q 2 )ApD + ifs(q 2 )AsD + 


ifD(q 2 )ADD 


(B.7) 


where D^/q 2 ) is the BW denominator defined in 
Eq. (IB.3I) and describing the fei(1235) :F shape, clds 
and <pos are the parameters describing the admix¬ 
ture of S and D waves in the amplitude of the 
6i(1235) =F decay and is the momentum fac¬ 

tor corresponding to the D wave in the &i(1235) =F 
decay. This factor can be defined for the intermedi¬ 
ate resonance with arbitrary integer spin J as 


rp(i450)(g 2 ) 


_m p (i450) kl ^ r p(1450) 

m p(1450) Po,ui % 


(B.6) 


where k, T is the same as in Eq. (IB.51) but computed in 
the ^(MSO)^ rest frame and fco, 7 r is calculated as k n 
but with \fq 2 = m p (i 45 o). The first term in Eq. (IB.61) 
corresponds to the ^(1450)^ —» decay while 

the second describes the —> ujtt p decay. We 

assume that the ^(1450)^ resonance decays to these 
final states with equal probabilities. For events with 
q 2 < (m u + TOtt) 2 , we use r p(1450) (g 2 ) = r p(1450) . 

We assume that the relative helicity phases <p± for 
the ^(ldSO)^ production are the same as for the off- 
shell p{ 770) T . This assumption does not contradict 
the common description of the matrix element be¬ 
cause of the validity of the factorization hypothesis. 
Since the typical values of q 2 are close to each other 
for the p(770) p and p(14A0) p , we can neglect the dif¬ 
ference between the appropriate FSI helicity phases 
for these resonances. 


Pj(P3tt ) 




(B.8) 


The form factors fs{q 2 ), fp{q 2 ) and /n(g 2 ) are de¬ 
termined in Appendix ICl and Ass, Aps, Ads, Asd, 
Apd and Add are the angular dependencies defined 
in Table IVII that correspond to the ui^ quantum 
numbers J p = 1 + . 

The parameters clds and epos are fixed at the 
values measured by the Brookhaven E852 collabo¬ 
ration j57||. There, the amplitude ratio was found to 
be |D/S| = 0.269 ± 0.013 and the phase difference 
4>ds = 0.18 ± 0.08 rad [57f. To relate the parame¬ 
ter clds to the ratio \D/S\, the helicity amplitude 
M ++ , corresponding to the positive helicities of the 
&i(1235) T and the to in the decay 6i(1235) =F —>• unr p , 
is written in terms of partial waves: 


M++ = M| + + M p + = S/s /.3 + D/s/6, (B.9) 

where M? + and M/? + are the terms corresponding to 
the S and D waves, respectively. To calculate these 
terms, we denote the polarization four-vectors of the 
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Table VI: Angular dependencies corresponding to the 
ojn^ quantum numbers J p = 1 + . L\ (L 2 ) is the relative 
orbital angular momentum between the D* ± and con^ 1 
(uj and 7r^~). The notations c a = cos a and s a = sin a 
are used. The angles 9, <j>, /3, ip, £ correspond to the urn 
angular basis. 


Li 

l 2 

^ l<2, 

S 

S 

CoC/3C£ SQC^CpS^ SQSSftS ifj 

+S0C ( /)SpG l jjC€ + cespc-^s^ 

P 

s 

CQSpS^S^ ~\~ SQS(j)SpCip SQC^SpS^C^ 

D 

s 

2CQCpC£ SQC(j)Cp S £ SQ S(f)Sp 

~\~SQC(j)S pCijjC^ 2 cq S pC-ijj S£ 

S 

D 

2cgCpC£ 2soC ( j ) CpS£ SqS^S pSij)~\~ 

~\~SffC(j)S pCijjC^ “l - CqS pCips^ 

P 

D 

2cq Sp S'ljjS^ ~\~ SQS<f)SpC-ijj SQC(f)sp S ip C£ 

D 

D 

4 CQCpc^ 2sQC(j)Cps £ sq s (j)S pSij}~\~ 

-j-sec^spCijjC^ — 2 ce sp c-ijj 


6i(1235) T and the w as e M and respectively. In 
such a case the terms are written as 


Since the typical values of q 2 for this resonant 
decay is close to the values corresponding to the 
p(770) T and p{ 1450) T amplitudes, we assume that 
the FSI helicity phases <j>± in this decay are the same 
as for the p(770) T and p(1450) :p contributions. The 
free parameters for this contribution are the relative 
magnitude and phase. 

B°(B°) -> Di(2430)°(Di(2430)°)w -» 

B°(B°) Di(2420)°(Di(2420)°)w 

The notations D[ and D\ for the Z?i(2430)° 
(Z?i(2430)°) and £>i(2420)° (14i(2420) 0 ), respec¬ 
tively, are used in this subsection. 

The observable D[ and Di states are not charge- 
conjugation eigenstates but rather the admixtures 
between the pure states with J p = ly 2 and Jf = 

where the quantum number j is the total an¬ 
gular momentum of the u quark Q. Mixing in the 
jj coupling scheme is written as 

Md i± = Daiq 2 ) ( ai / 2e ^ V2 sinujM i/2± + 
a 3/2 e% ^ 3/2 cos we _iV M 3/2± ^ , 

M D ' 1± = D ^ q 2^ ( a l/ 2 e ^ 1/2 COSUlM 1/2 ± - 

a 3/2 el<l>3/2 smive iV> M 3 / 2 ±\ , (B.13) 


M+ + = m 2 bl e+v + * M , 

Mf + = -a DS e^ DS Pi(po,u 1 )£+v+^,( B.10) 

where p 2 = m 2 in Pi(po,u>)- Taking into account 
that e+v + * A ‘ = —1, we have 


ods 


1 \D\ 

V2 p i(Po,u:) |S| ’ 


(B.ll) 


and obtain clds = 5.2 ± 0.3. 

The width r bl ( q 2 ) for events with q 2 > (m^-t-m^) 2 
is parametrized via the 6i(1235) =F —> unr T decay: 


r 6 l (g 2 ) 


mbj_Pu_ 

^Po,. &1 

nib^siPu) + ^dsPHp^BKp^) 

m t + 2a DS P l(P0,u) 

(B.12) 


where ro bl(1235) = (1230 ± 3) MeV/c 2 , T bl = (142 ± 
9) MeV/c 2 !3j and the factor 2 accounts for the nor¬ 
malization of the D wave relative to the S wave. For 
events with q 2 < (jn u +m T ) , we use r bl (< 7 2 ) = F bl . 


where u> and <p are the mixing angles, 0 , 1 / 2 , a 3 / 2 , 
4 > 1/2 and <f> 3/2 are the relative magnitudes and phases 
between the pure matrix elements M\ / 2 ± and M 3 / 2 ±, 
which correspond to the Jj’ = ly 2 and = 1^, 2 
quantum numbers, respectively. 

The pure matrix elements Mi/ 2 ± and M 3 / 2: t are 

Mi/ 2± = im 2 D , i Bs(p D -)[±fp(q 2 )Aps + 

ifs{q 2 )-Ass + */ z >(<? 2 )- 4 ds ], 
a 'h/ 2 ± = —i P l(PD*)BD{pD*)[±fp(q 2 )ApD + 

ifs{q 2 )-AsD+ifD(q 2 )ADD], (B.14) 

where P\(pd*) is defined in Eq. (1B.8|) with p 2 = 
m 2 D , and the angular dependencies Mss, Aps, Ads, 
Asd 1 Apd and Add have the form shown in Ta¬ 
ble IVII except for one feature: the angular basis (0, 
</>, /3, V’j 0 describes here the D** production. The 
transition form factors fs(q 2 ), fp(q 2 ) and /p(g 2 ) 
are given in Appendix [Cl Since the mixing effect 
is predicted and confirmed to be small, we use in 
Eq. (IB.14[) the physical mass mp>' instead of the mass 
of the pure j = 1/2 state. 
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The g 2 -dependent widths Tp 1 (q 2 * ) and T D[(q 2 ) of 
the D i and D[ states are parametrized via their de¬ 
cays to D 


r D ,(q 2 ) = -^B 2 {p D ,) 

r -Di(? 2 ) = —F= B d(PD*) 


D’ 


Pq,d* 

p?(pd* 


Pd* 




P{ {PO,D*)po,D> 


r Du 
(B.15) 


where P\{pd*) and Pi(po,d*) are defined as in 
Eq. (fTTTni . m Dl = (2421 ± l)MeV/c 2 , T Dl = 
(27 ± 3) MeV/c 2 , m D[ = (2427 ± 36)MeV/c 2 and 
V D[ = (384 ± 117) MeV/c 2 are fixed Q. 

The free parameters describing these resonant am¬ 
plitudes are the mixing angles w and p, the relative 
magnitudes and phases aq/ 2 , a 3 / 2 , 4 >i /2 and <(> 3 / 2 , and 
the relative normalizations and phases defined in the 
/s(g 2 ), fp(q 2 ) and fp(q 2 ) form factors. 

B°(B°) -» D 2 (2460)°(D 2 (2460)°)u; -> 

The notation D 2 for the £>2(2460)° ( 742 ( 2460 )°) is 
used in this subsection. 

Since the D | state has the quantum numbers 
J p = 2 + , the P and F waves describing its pro¬ 
duction violate P- and (7-parities, and thus have the 
additional phase 7 r /2 in comparison with the D wave. 
The resonant matrix element is 


M d *± = 


Pd* b d{pd*] 


fD(q 2 )App ± 


D D *(q 2 ) 
ifp(q 2 )Ap D ± ifp(q 2 )AFD 


(B.16) 


where Apd, Add and Afd are the angular depen¬ 
dencies describing each partial wave and shown in 
Table fVTTl and the transition form factors fp{q 2 ), 
fp(q 2 ) and /^(g 2 ) are parametrized in Appendix [Cl 
The g 2 -dependent width Td* (g 2 ) is determined via 
decays of the H 2 to D and D ± n T with the 
probabilities of 40% and 60%: 


r D *(q 2 ) 


2 y /7 P 5 p . 

5 rn D * p 5 0 D , 


Bd(pd*)Fd* + 


3 m Dj kp 
5 y/g 2 D 


Bp(kp)r D *, 


(B.17) 


where kp is the D ^-meson momentum magnitude 
in the D 2 —> D ± tt^ decay computed in the D 2 rest 
frame, ko^p is the same momentum when y/g 2 = 
m D * = (2463 ± 1) MeV/c 2 , and T p* = (49 ± 
1) MeV/c 2 3. 

The free parameters, describing this tensor contri¬ 
bution, are an overall magnitude and phase as well 
as normalizations and relative phases of the partial 
wave form factors in the matrix element. 


Table VII: Angular dependencies corresponding to the 
D f± quantum number J p = 2 + . L\ (L 2 ) is the 
relative orbital angular momentum between the D* ± 7 r =F 
and oj ( D* ± and 7 ^). The notations c a = cos a and 


s a = sin a 
to the D** 

are used. The angles 9, cp, /?, ip, £ correspond 
angular basis. 

Li 

l 2 

vALi 1/2 

P 

D 

Cq S S'ljj S2£ + SQC'ftSpS^Z S q S (fi S ft C-ijj 

D 

D 

SQSifiS 0 Sijj S0 C(f)S/3Cijj C£ 

F 

D 

3 / 200 SfiSijjS2£ SqC(J)S/ 3SijjC2£ SqS 


APPENDIX C: PARTIAL WAVE FORM 
FACTORS 


In this section, we obtain full expressions of the 
partial wave form factors used in the resonant matrix 
elements. The symbols and Pd*,b are used in 
this section for the magnitudes of the three-momenta 
of the w decay product system and D * ± in the B 
meson rest frame, respectively. 

The partial wave form factors fs{q 2 ), fp(q 2 ) and 
fp(q 2 ) describing the ujtt t resonance production in 
Eqs. (IB.21) and (IB. 71) can be expressed in terms of 
three helicity amplitudes [H 0 (q 2 ) and H±{q 2 )\, which 
correspond to three polarization states of the D * ± 
(one longitudinal and two transverse), and two trans¬ 
verse helicity phases <j>± defined relative to the lon¬ 
gitudinal amplitude H 0 (q 2 ): 


fs(q 2 ) = 
fp(q 2 ) = 
fo{q 2 ) = 


q 2 H+(q 2 )e^+ + i7_(g 2 )e^- + H 0 (q 2 ) 

3 

g2"i7+(g 2 )e^+ - i7_(g 2 )e^- 
2 ’ 
^H+{q 2 )e^+ + H.{q 2 )e^~ - 2 H 0 {q 2 ) 


V6 


(C.l) 


Here, the additional factors y/g 2 /2, y/g 2 /3 and 
Vf/G are introduced to take into account the g 2 - 
dependent vertex of the W 7 r production in the fac¬ 
torization assumption and the relative normalization 
fractions of the angular dependencies shown in Ta¬ 
bles 0 and ED The helicity amplitudes H 0 (q 2 ) and 
H±(q 2 ) can be written in terms of three invariant 
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form factors A\(q 2 ), A 2 {q 2 ) and V(q 2 ) lllil: 


p 2 = 0.79 ±0.06. (C.6) 


H 0 (q 2 ) 


H±{q 2 ) 


(m|> - m 2 D , - q 2 )(m B + tod*) 


2 to d * V 7 ? 2 


2p 


2 

D*,B 


m. 


m D * \[q 2 (m B ± m D *) 

- (m B + m B *)Ai(q 2 ) ± 


A 2 (q 2 


Ai(q 2 


2pD*,Bm B 2 

—± - v(q )■ 

m B ± mu* 


(C.2) 


The invariant form factors Ai(q 2 ), A 2 {q 2 ) and V(q 2 ) 
describe the B°(B°) —»■ D * ± transition and can be 
related to the Isgur-Wise function h{w) ( w is the in¬ 
variant four-velocity transfer) under the assumption 
of heavy quark symmetry 


A^q 2 ) 

A 2 (q 2 ) 

V(q 2 ) 


1 - 


r 2 

Ri 


{m B ± m D *) 2 
m B ± mp 


y/m B m D * 
m B ± tod* 
y/m B m D * 


h(w), 

h(w), 


m B ± tod* 
2 y/m B m D * 


h(w), 


(C.3) 


where R\ and R 2 are the relative factors, and the 
Isgur-Wise function h(w) can be parametrized as [1, 


h(w) = l-p 2 (w-l) (C.4) 


with 


to! + to!, — q 2 

2 TObTOu* 


(C.5) 


The values R\ : R 2 and p 2 used in our analysis were 
measured by the BABAR collaboration in the B° —> 
D* + e~ u e decay 471: 


i?i = 1.40 ±0.06, 
R 2 = 0.87 ±0.04, 


The partial wave form factors fs(q 2 ), fp{q 2 ), 
fp{q 2 ) and /f(? 2 ) describing the D** resonance pro¬ 
duction and introduced in Eqs. (IB. 141) and (IB.161) 
contain the momentum dependencies, corresponding 
to the definite angular orbital momenta of the decay 
products in the B meson rest frame. These depen¬ 
dencies can be explicitly extracted. The form fac¬ 
tors for the ±>i(2430)° (£±(2430)°) and ±>i(2420) 0 
(±>i(2420)°) production can be written as [36|] 


fs(q 2 ) = --j=m 2 B B s {q 2 )e l(,,s , 
fp{q 2 ) = ^m B p 3 iv, B Bp(q 2 )e lcl>F , 

f D (q 2 ) = ^ Pi(p3n,D**)B D (q 2 ), (C.7) 

where Rs and Rp (<ps and <j>p) are magnitudes 
(phases) of S- and P -wave amplitudes defined rel¬ 
ative to D -wave amplitude, ±i(p 37 r,D**) is defined in 
Eq. (IB. 81) and p 37r , £>»* is the three-momentum mag¬ 
nitude of the w decay product defined in the D* ± Tr T 
rest frame. Similar expressions can be written for 
the D2(2460) 0 (-£>2(2460)°) production [36j : 


fp(q 2 ) 
f D {q 2 ) 
fp(q 2 ) 


R P rn‘ B p3 7T:B 
a/3 to d * 


Bp(q 2 y* p , 


Rp m BP 2 3„, B 

a/2 tod* 


Bp(q 2 y 


‘D 


P 2 (P3tv,D**)P3tt,B 

2 tod* 


B F (q 2 ), 


(C.8) 


where Rp and R B {<fip and 4>d) are magnitudes 
(phases) of P- and D -wave amplitudes defined rel¬ 
ative to F -wave amplitude, P 2 (j> 3 tt,d**) is defined in 
Eq. (IB. 81) with J = 2 and tod* is the mass of the 
±>2(2460)° (±>|(2460)°) state. 
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